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Abstract 

We present an improved measurement of the branching fraction for the electroweak penguin 
process B — > Xs£'^i~, where i is an electron or a muon and Xs is a hadronic system containing 
an s-quark. The measurement is based on a sample of 152 x 10^ T(4S') ^ BB events collected 
with the Belle detector at the KEKB e~^e~ asymmetric-energy collider. The Xg hadronic system 
is reconstructed from one or Kg and up to four pions, where at most one pion can be neutral. 
Summing over both lepton flavors, the inclusive branching fraction is measured to be B{B 
Xs£+e-) = (4.11 ±0.83(stat)l[J:^^(syst)) x 10"^ for m(£+r) > 0.2 GeV/cP. 

PACS numbers: 13.20.He, 12.15.Ji, 14.65.Fy, 14.40.Nd 
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INTRODUCTION 



In the Standard Model (SM), the rare decay B Xsi~^i~ proceeds through a 6 ^ s£~^£~ 
transition, which is forbidden at tree leveL Such a flavor-changing neutral current (FCNC) 
process can occur at higher order via electroweak penguin and W^W~ box diagrams. The 
b si~^i~ transition therefore allows deeper insight into the effective Hamiltonian that 
describes FCNC processes and is sensitive to the effects of non-SM physics that may enter 
these loops; see, for example, Refs. 

Recent SM calculations of the inclusive B Xsi~^i~ branching fractions predict 13{B — > 



XsC+e-) = (6.9 ± 1.0) X 10"' 
Xsfi+fi-) = (4.2±0.7)xl0-6 
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4.2 ± 0.7) X 10-6 for m{£+£~) > 0.2 GeV/c^] and B{B 



Both the Bel 



exclusive B Kl^l" and K*£+£~ decays 
Xs£+£- decay QQ. 



e and BaBar Collaborations have observed 
3 a.d W .ea.u.d .elusive B 



U-uau 



In this analysis, we study the inclusive B Xs£^£~ process by semi-inclusively recon- 
structing the final state from a pair of electrons or muons and a hadronic system consisting 
of one or and up to four pions, where at most one pion can be neutral. This 
semi-inclusive-reconstruction approach jlO|] allows approximately 53% of the full inclusive 
rate to be reconstructed. If the fraction of modes containing a is taken to be equal to 
that containing a Kg, the missing states that remain unaccounted for represent ~30% of 
the total rate. We require the hadronic mass for the selected final states to be less than 
2.0 GeV/c^ to reduce combinatorial background. We correct for the missing modes and 
the effect of the hadronic mass cut to extract the inclusive B Xs£^£~ decay rate for 
m{£+£-) > 0.2 GeV/c2. 

The measurement of inclusive B —>■ Xs£'^£~ decay in this paper updates and supersedes 
our previous publicationP| described above, which was based on a sample of 65 x 10^ BB 
pairs. The measurement reported here is currently the most precise. 

THE BELLE DETECTOR AND DATA SAMPLE 

We use a data sample collected on the T(4S') resonance with the Belle detector at the 



KEKB e~^e asymmetric-energy collider (3.5 GeV on 8 GeV) 



This sample comprises 



152 X 10^ B meson pairs, corresponding to an integrated luminosity of 140. Ofb ^. A detailed 
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description of the Belle detector can be found elsewhere A three-layer silicon vertex de- 
tector (SVD) and a 50-layer central drift chamber (CDC) are used for tracking and particle 
identification for charged particles. An array of aerogel threshold Cerenkov counters (ACC) 
and time-of-fiight scintillation counters (TOF) are used for the charged particle identifica- 
tion. An electromagnetic calorimeter comprised of Tl-doped Csl crystals (ECL) measures 
the energy of electromagnetic particles and is also used for electron identification. These 
detectors are located inside a superconducting solenoid coil that provides a 1.5 T magnetic 
field. An iron flux-return located outside of the coil is instrumented with resistive plate 
counters to identify muons (KLM). 

Particle identification for the final state particles e^, fi^, K^, K^, vr^ and 7r° is important 
for this analysis. Electron identification is based on the ratio of the cluster energy to the 
track momentum (E/p), the specific energy-loss measurement (dE/dx) with the CDC, the 
position and shower shape of the cluster in the ECL and the response from the ACC. Muon 
identification is based on the hit positions and the depth of penetration into the ECL and 
KLM. Electrons and muons are required to have lab-frame momenta greater than 0.4 GeV/c 
and 0.8 GeV/c, respectively. For the muon identification, we also apply a kaon veto to select 
good muon candidates. Bremsstrahlung photons from electrons are recovered by combin- 
ing an electron with photons within a small angular region around the electron direction. 
Charged kaon candidates are selected by using information from the ACC, TOF and CDC. 
The kaon selection efficiency is 90% with a pion to kaon mis-identification probability of 6%. 
After selecting the electron, muon and charged kaon candidate tracks, the remaining charged 
particles are assumed to be charged pions. candidates are reconstructed from pairs of 
oppositely-charged tracks with \m{TT^7T^) —m{Kg ) \ < 15 MeV/c^. We impose additional Kg 
selection criteria based on the distance and the direction of the Kg vertex and the impact 
parameters of the daughter tracks. We require the charged tracks except for those used in 
the Kg reconstruction to have impact parameters with respect to the nominal interaction 
point of less than 1.0 cm in the radial direction and 5.0 cm along the beam direction. Neutral 
pions are required to have lab-frame energy greater than 400 MeV, photon daughter energies 
greater than 50 MeV, and a 77 invariant mass that satisfies |m(77) — m(7r*^)| < 10 MeV/c^. 
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ANALYSIS OVERVIEW 



In this analysis, we reconstruct the inclusive B Xsi^i~ decays with a semi-inclusive- 
reconstruction technique from a pair of electrons or muons and one of 18 reconstructed 
hadronic states. Here the hadronic system consists of one or and up to four pions 
(at most one pion can be neutral). Compared to a fully inclusive approach, this method 
has the advantage of having strong kinematical discrimination against background by using 



the beam-energy constrained mass Mbc = y E'^eam " p'b energy difference AE = 

Eb — Ebeam, whcrc Ebeam IS the beam energy and Eb (pb) is the reconstructed B meson 
energy (3-momentum). All quantities are evaluated in the e~^e~ center-of-mass system (CM). 

In addition to the discrimination, further background suppression to reduce the large com- 
binatorial backgrounds is necessary. The main contribution to the combinatorial background 
comes from semileptonic decays in BB events. In these events, B — > Xgi^i^' candidates 
are reconstructed with the decay products from both BB mesons. This background has 
a significant amount of missing energy due to the neutrinos from the semileptonic decays. 
Another contribution to the combinatorial background comes from continuum events, which 
are effectively suppressed with event-shape variables. 

There are two background sources that can peak in Mbc and AE. The first comes from 
B J/ipX and B ip{2S)X decays with J/ip{ip{2S)) i'^i^. This peaking background 
is efficiently removed with cuts on the dilepton mass m{i~^i~). The resulting veto sample 
provides a large control sample of decays with a signature identical to that of the signal. The 
second comes from B K^{Kg)mT {n > 1) decays with misidentification of two charged 
pions as leptons. We estimate these peaking background contaminations, then subtract 
them from the signal yield. 

For the B Xsi~^£~ event simulation, we use EVTGEN for the event generator, 
JETSET to hadronize the system consisting of a strange quark and a spectator quark, 
and GEANT for the detector simulation. In the event generation, B Xsi~^i~ events 
are produced with a combination of exclusive and inclusive models. In the hadronic mass 
region of miXs) < 1.1 GeV/c2, exclusive B K^*H+t decays are generated according 
to Refs. 0, Q], where the relevant form factors are computed using light-cone QCD sum 
rules. In the region m{Xs) > 1.1 GeV/c^, event generation is based on a non-resonant model 
following Refs. ^, 12 1 and the Fermi motion model of Ref. jlS^ . 
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EVENT SELECTION 



Events are required to have a well determined primary vertex, be tagged as multi- 
hadron events, and contain two electrons (muons) having lab-frame momenta greater than 
0.4 GeV/c for electrons and 0.8 GeV/c for muons. Dilepton candidates are selected for e^e~ 
or fi^fi~ pairs. Both leptons are required to originate from the same vertex and satisfy the 
requirement \Az\ < 0.015 cm. Here, Az is the distance between the two leptons along 
the beam direction; the ^-coordinate of each lepton is determined at the point of closest 
approach to the beam axis. 

Charmonium backgrounds are reduced by removing B candidates with a dilepton mass in 
the ranges -0.40 GeV/c^ < M^e(y) - Mj/^ < 0.15 GeV/c^, -0.25 GeV/c^ < M^^ - Mj/^ < 
0.10 GeV/c2, -0.25 GeV/c^ < Mee{^)- M^i2S) < 0.10 GeV/c^, and -0.15 GeV/c^ < M^^- 
M^{2S) < 0.10 GeV/c^. If one of the electrons from a J/iIj or '?/'(2S') decay erroneously 
picks up a random photon in the Bremsstrahlung-recovery process, the dilepton mass can 
increase sufficiently to evade the above cuts. Therefore the charmonium veto is applied to the 
dilepton mass before and after Bremsstrahlung recovery. Using the simulation, we estimate 
the remaining peaking charmonium background to be 1.20 ± 0.28 events and 1.33 ± 0.21 
events for e^e~ modes and ^~ modes, respectively. 

The potential peaking background from B decays, followed by conversion of the 

photon into an e^e~ pair in the detector material, and 7r° Dalitz decay is a concern for the 
e^e~ modes only. We remove this background by requiring m{e~^e~) > 0.2 GeV/c^. 

Using the pair, B Xsi~^i~ candidates are formed by adding either a K"^ 

or a and up to four pions, but no more than one 7r°. In this manner, eighteen 
different hadronic topologies are considered: K^, K^n^, K^tt^, i^^Tr^Tr*^, K'^tt^tc^, 

K^TT^n^TT^, K'^TT^n^TT^, K'^TT^TT^TT^TT^, K ^ Tl^ TV^ TT^ TT^ , K^, K^-n^ , K'^gTr^ , K'^TT^TT^, 

K'^sTT^n^, K'^gn^n^TT^, K'^gn^n^n^, K'^sn^n^n^n^, and K'^sTi^n^n^n^. 

After forming the B —>■ Xgi'^i" candidates, we carry out the background suppression. 
The largest background sources are random combinations of dileptons with a kaon and 
pions that originate from continuum qq {q = u, d, s, c) production or from semileptonic 
B decays. We reject the qq background with a Fisher discriminant |3| (^sfw) based 
on a modified set of Fox- Wolfram moments [2^ that differentiate the event topology, by 
applying a cut J^^sfw > —1-0. In the semileptonic B decay background, both B mesons 
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decay into leptons or two leptons are produced from the b c s,d decay chain. We 
combine the missing mass (Mmiss) and the total visible energy (-Evis) into another Fisher 
discriminant (J-'miss) to reject the B decay background. The missing mass is defined as 
Mmiss = \/{'^Ebeam - Ei)"^ - | EPiP, whcre Ebeam IS the CM beam energy and Ei (pi) is 
the reconstructed energy (3-momentum) for the charged particle or the photon in the CM 
frame. The sum runs over all the charged particles with a pion mass hypothesis and all 
photons. The total visible energy is also calculated using the charged particles and photons. 
We require the selection cut of J?-miss > 1-2. 

We further reduce the background using AE and XBvtx with the selection criteria: 
-0.10 GeV/c < AE < 0.05 GeV/c for the electron mode (-0.05 GeV/c < AE < 
0.05 GeV/c for the muon mode), and x\vtx/^E)F < 10.0. Here, x%vtx is the fitted chi 
squared of the B vertex constructed from the charged daughter particles except for the 
daughter tracks. We reject candidates with the invariant mass greater than 2.0 GeV/c^. 
This condition removes a large fraction of the combinatorial background while retaining 99% 
of the signal. 

At this stage, there is an average of 1.6 B candidates per event in the signal simulation. 
We retain only the B candidate with the largest signal likelihood. We select the following six 
variables as the background-discrimination variables, and calculate the likelihood functions 
based on the distributions of the variables: AE, AE^'^^, x%vtx^ cos 9b, ^sfw, and J-miss- The 
energy difference AE^^^ is formed by combining all charged tracks and neutral calorimeter 
clusters not included in the B candidate, and cos 9b is the cosine of the B flight direction 
with respect to the e~ beam direction in the CM frame. 

The signal probability density functions (PDFs) are determined by applying flts to each 
distribution for signal MC events. For AE and x'Bvtx^ the real charmonium- veto- 

event distributions to determine the PDFs, because we observe some discrepancy between 
the signal MC and the real charmonium-veto events. In the charmonium-veto-event distri- 
butions, we subtract the background shapes obtained by the Mbc side-band region. The 
normalization of the subtracted background events is determined by the number of back- 
ground events in the Mbc signal region, which is estimated by fltting to the Mbc distribution 
with a Gaussian (signal) and an ARGUS (BG) function. For AE, we use distinct PDFs 
for the electron and muon modes. For the background PDFs, we use background MC events. 

The variables AE, AE^^^ , and J^mi^s are effective at rejecting BB background, especially 
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XsJ/\\f data 
Xsll MC 
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Xseia data 
Xseia MC 



^ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I ri I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.1 02 0.3 04 0.5 0.6 0.7 OS 0.9 1 0.1 02 0.3 04 0.5 0.6 0.7 OB 0.9 1 



Likelihood Ratio 



Li[<eliliood Ratio 



FIG. 1: The likelihood ratios for (a) signal and (b) background events. The red (blue) and black 
histograms in the signal (background) likelihood distributions show signal MC (MC B — > XgC^^^ 
events) and the real charmonium veto sample (real B XsC^fi^ events), respectively. 



for events with two semileptonic decays, which have large missing energy. For continuum 
suppression, the event-shape variable J-sfw and cos 9b are useful. XBvtx is effective to reject 
the random combinatorial background in the high multiplicity modes. 

We then calculate the likelihoods Cs,b = Y[^=iPs,b where are the PDFs for the 
background-discrimination variable i for the signal (S) and the background (B), respectively. 
Then we obtain the final discriminating variable, the likelihood ratio TZ = Cs/{Cs + Cq). 
Only the B candidate with the largest TZ value is retained. We find that the correct candidate 
is reconstructed in 84% of the events. 

In order to check the obtained likelihood ratio, we compare the real charmonium-veto 
events and signal MC for the signal likelihood ratio, and the real and background MC 
B XgC^fi^ events that are selected using the nominal selection criteria but requiring that 
the two leptons have different flavor, for the background likelihood ratio. Figures ^a) and 
(b) show the likelihood ratios for the signal and background events, respectively. We observe 
good agreement in both flgures. 

The flnal suppression of the combinatorial background is achieved with the likelihood ratio 
TZ. Using the simulation, the cut on TZ is optimized to maximize the statistical signiflcance of 
the signal. This optimization is performed in the two m{Xg) regions m{Xs) < 1.1 GeV/c^, 
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and 1.1 GeV/c^ < m(Xs) < 2.0 GeV/c^, resulting in the cuts 7^ > 0.3 and 0.9. 

After applying all selection criteria, a sample of 155 B — > Xge'^e" and 112 B — > 
candidates remains in the signal Mbc region. Here we define the signal M^c region as 
5.27 GeV/c^ < M^c < 5.29 GeV/c^. According to the simulation, the background remaining 
at this stage of the analysis consists mostly of BB events (80% and 72% of the total back- 
ground in the electron and muon channels, respectively). Using the signal MG simulation, 
the probability to select the correctly reconstructed candidate is estimated to be 91%. 

MAXIMUM LIKELIHOOD FIT 

We perform an extended, unbinned maximum likelihood fit to the Mbc distribution in the 
region Mbc > 5.2 GcV/c^ to extract the signal yield as well as the shape and yield of the 
combinatorial background. The likelihood function C is expressed as: 

^ = n[(^-. + N,eak)rr + E NBovr] 

■ i=l BG 

Ntotal-BG = Npc + Ncf + Ncomb 

E NBGVf'' = N,cVr + NcfVf + Nco^bVr"' 

BG 

where N and i denote the total number and index of candidate events, respectively. Ngig, 
Npeak, Npc, Ncf, and Ncomb represent the yields of the signal, peaking background, combina- 
torial background from peaking background, cross-feed events from the mis-reconstructed 
B Xgi'^i" decays, and combinatorial background, respectively. V-^^ is the signal PDF. 
We use the same PDF for signal and peaking background events. Pf, Vf, and ^ are 
the background component PDFs for combinatorial background from peaking background, 
cross feed, and combinatorial background, respectively. 

The signal PDF Vp^ is described by a Gaussian for /x^/i" modes as well as for e+e" 
modes, since the Bremsstrahlung recovery and selection procedure for e"*"e~ modes lead to a 
negligible radiative tail in the Mbc distribution. The Gaussian shape parameters, the mean 
and the resolution a values are determined from fits to the sum of a Gaussian and an ARGUS 
function for the real charmonium-vcto data sample. The fits results in a signal Mbc mean 
and resolution, respectively, of m^^g = 5279.31±0.05 MeV/c^ and a^ig = 2.62±0.04 MeV/c^ 
for the e+e" modes, and rUsig = 5279.03 ± 0.04 MeV/c^ and a^ig = 2.53 ± 0.04 MeV/c^ for 
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the jjL^jjL~ modes. In the simulation, the Gaussian fit results for the M^c distributions for 
correctly reconstructed signal are in agreement with the shape parameters extracted from 
the fits to the charmonium-veto sample. The signal yield Ng^g is a free parameter in the 
likelihood fit. 

The charmonium peaking background is estimated from the simulation to be 1.20 ± 0.28 
events in the electron modes and 1.33 ± 0.21 events in the muon modes. The charmonium 
peaking background PDF is the same as that for signal since the signal PDF is extracted 
from the charmonium-veto data sample. 

The size and shape of the hadronic peaking BB background component arising from 
B D^*^n7c [n > 0) decays with misidentification of two charged pions as leptons are derived 
directly from the real data by performing the analysis without the lepton identification 
requirements. By fitting the M^c distribution to the sum of Gaussian and ARGUS functions, 
we get the mean and the resolution of nipeak = 5279.16 ± 0.04 MeV/c^ and <Jpeak = 2.60 ± 
0.02 MeV/c^. Taking the n to i misidentification rates into account, the remaining hadronic 
peaking background is estimated to be Nh-peak = 0.03±0.001 events for the e~^e~ modes and 
Nh-peak = 1-78 ± 0.05 events for the fi'^fi^ modes. Here we use the momentum and polar- 
angle dependent misidentification rate. The average misidentification rates for electrons 
and muons are 0.08% and 0.92%, respectively. In the likelihood fit, we also use the same 
signal PDF for the hadronic peaking-background PDF, because the fitted Gaussian shape- 
parameter values are consistent with those of the charmonium-veto data sample, and Np^ak 
is fixed to the estimated values. Figures El^a) and (b) show the Mbc distributions for the 
MC charmonium events for the electron and muon modes, and (c) and (d) show the Mbc 
distributions for real B Xgh^h" events for the electron and muon modes, respectively. 

The background component PDFs, Pf^, V^^ , and 'P^°'^'^ are given by an ARGUS shape. 
They describe the combinatorial contribution from peaking background events, from cross 
feed events, and from continuum events and BB events, respectively. The ARGUS cutoff is 
determined by the beam energy in the T(4S') rest frame, Ebeam = 5.290 GeV. The values 
of the ARGUS shape parameter for each background component are determined from the 
peaking background M^c distribution shown in FigI21 {Vf^) , from incorrectly reconstructed 
signal MC events (V^-^), and from the fit to the real B Xge^ir^ events selected using the 
nominal selection criteria but requiring that the two leptons have different flavor (Pf 
We fix the three ARGUS shape parameters, Np^, and N^f. The yield Ncomb is taken as a 
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(a)Xs charmonium : e+e- (b) Xs charmonium : 




Mbc (GeV/c2) Mbc (GeV/c2) 



FIG. 2: Mbc distributions for MC charmonium events for (a) electron and (b) muon modes, and 
real B Xsh^h~ events for (c) electron and (d) muon modes. 

free parameter in the hkehhood fit. 
RESULTS 

Using the fit parameterizations described above, we fit the Mbc distributions for the 
selected B X^e^e^ and B —>■ Xsfi~^fi~ candidates separately and obtain the results 
shown in Figure El The fit results are summarized in Table |l] The statistical significance 
is S = 2{\n Cmax — ln£°ja2,), where Cmax represents the maximum likelihood for the fit 
and C^^^ denotes the maximum likelihood for a different fit when the signal yield is fixed 
at Nsig = 0. The B — * Xsi~^i~ signal yield presented in Table H] is the sum of the B — »• 
Xse'^e~ and B Xg^^ ^~ signal yields. A separate fit to the combined electron and muon 
channels gives a comparable result. FigureEfd) shows the M^c distribution for B XsC^jj^ 
candidates. Applying the ARGUS fit to the M^c distribution, there is no evidence for a 
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peaking background as expected. 

Figures|3fa) and (b) show the distributions for the hadronic mass Mx^ and = Mj+^- for 
the electron and muon channels combined, but they are obtained by performing the nominal 
likelihood fit in separate Mx^ and regions. Figure El^a) indicates that the observed signal 
includes contributions from final states across a range of hadronic masses, including hadronic 
systems with a mass above that of the fC*(892). 

The branching fraction B for the signal is calculated from 

■ 

where N^j^ = (152.0 ± 0.7) x 10^ is the number of BB pairs produced in 140.0 fb~^ and e 
is the signal efficiency. 



SYSTEMATIC UNCERTAINTIES 



Systematic uncertainties are of two different types: those that affect the extraction of the 
number of signal events and those that affect the calculation of the branching fraction. The 
systematic uncertainties are summarized in Table IHI 

Uncertainties affecting the extraction of the signal yield are evaluated by varying the 
signal Gaussian parameters (mean and width) and the background shape parameter within 
±lcr of the measured values from the charmonium veto data (signal) and the real B — »• 
XsC^/x^ events (background). 

We estimate the uncertainties in the peaking-background shape and cross-feed events 
by comparing the signal yields obtained with and without the corresponding PDFs in the 
unbinned maximum likelihood fit to M^c- 

TABLE I: Results of the fit to the data: number of the signal candidates in the signal box, obtained 
signal yield, peaking backgrounds (fixed in the fit), and statistical significance. 



Mode Candidates Ngig ^peak Signif. 

Xse+e- 155 31.8 ± 10.2 1.24 ±0.28 3.6 

Xsfi-^fi- 112 36.3 ± 9.3 3.11 ±0.22 4.7 

Xsi+i~ 267 68.4 ± 13.8 4.35 ± 0.36 5.8 
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TABLE II: Summary of fractional systematic uncertainties (in %). The uncertainties in extract- 
ing the signal are presented first and those related to the signal efficiency and BB counting are 
presented second. 



Source Xg 


e+e" Xs 




Signal shape 


±1.4 


±0.5 


BG shape 


±7.8 


±4.7 


Peaking background statistics 


±0.9 


±0.6 


Peaking background shape 


±4.3 


±2.1 


Cross-feed events 


±4.1 


±2.2 


Signal yield total 


±9.9 


±5.7 


Tracking efficiency 


±3.5 


±3.5 


Lepton identification efficiency 


±1.0 


±2.4 


Kaon identification efficiency 


±0.8 


±0.8 


TT^ identification efficiency 


±0.6 


±0.5 


Kg efficiency 


±0.7 


±0.8 


vr*^ efficiency 


±0.3 


±0.3 


TZ cut efficiency 


±5.4 


±1.5 


Detector model subtotal 


±6.6 


±4.7 


Fermi motion model 


+6.5 
-2.4 


+6.1 
-2.3 


B{B K£+e-) 


±6.4 


±6.8 


B{B K*i+i-) 


±7.0 


±7.8 


K*-Xs transition 


±4.5 


±4.7 


Hadronization 


±8.5 


±8.2 


Missing modes 


±4.5 


±4.4 


Signal model subtotal 


+15.6 
-14.4 


+15.9 
-14.9 


Monte Carlo statistics 


±1.6 


±1.5 


BB counting 


±0.5 


±0.5 


Efficiency and N^^ total 


+17.0 
-15.9 


+16.6 
-15.7 


Total 


+19.7 
-18.8 


+17.6 
-16.7 
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FIG. 3: Distributions of M^^ for selected (a) B X,e+e", (b) B Xsfi^fi', (c) B XjM' 
{i = e,/i), and (d) B X^e^fr^ candidates. The red lines represent the result of the fits, and the 
green and blue lines represent the peaking and combinatorial background components under the 
signal peaks, respectively. 



Uncertainties affecting the signal efficiency originate from the detector modeling, from the 
simulation of signal decays, and from the estimate of the number of B mesons in the sample. 
By far the largest component is that due to the simulation of signal decays, discussed in 
detail below. 

The detector modeling uncertainty is sensitive to the following uncertainties determined 
from the data: the uncertainty in the tracking efficiency of 1.0% per track; the uncertainty in 
the charged-particle identification efficiency of 0.5% per electron, 1.2% per muon, 1.0% per 
kaon and 0.8% per pion; and the uncertainty in the reconstruction efficiency of 4.5% per Kg 
and 3.3% per tt^. The efficiency of the likelihood ratio cut, which suppresses combinatorial 
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FIG. 4: Distributions of number of signal events as a function of (a) hadronic mass Mx^ and (b) 
q'^ = for electron and muon channels combined for data (points) and Monte Carlo signal 

(histogram). The vertical error bars represent statistical errors only. 



background, is checked with the charmonium-veto sample and the level of discrepancy with 
the simulation is taken as the corresponding uncertainty. 

The dominant source of uncertainty arises from modeling the signal decays. Parameters 
of the Fermi motion model are varied in accordance with measurements of hadronic moments 
in semileptonic B decays and the photon spectrum in inclusive B —>■ Xsj decays . 
The fractions of exclusive B Ki~^i~ and B K*i'^£~ decays are varied according to 
experimental 0, Q] and theoretical uncertainties respectively. The transition point in 
m{Xs) between pure K*[^[~ and non-resonant Xgl^l' final states is varied by ±0.1 GeV/c^. 

The non-resonant Monte Carlo event generator relies on JETSET to fragment and 
hadronize the system consisting of a final state s quark and a spectator quark from the 
B meson. Since the signal efficiencies depend strongly on the particle content of the fi- 
nal state, uncertainties in the number of charged and neutral pions and in the number of 
charged and neutral kaons translate into a significant uncertainty in the signal efficiency (for 
m{Xs) > 1.1 GeV/c^). 

The ratio between the generator yield for decay modes containing a Kg and that for modes 
containing a charged kaon is varied according to 0.50 ± 0.11, to allow for isospin violation in 
the decay chain. The ratio between the generator yield for decay modes containing one 7r° 
meson and that for modes containing none is varied according to 1.0 ±0.22. Uncertainties in 
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the ratios are set by the level of discrepancy between B ^ J / ipX real data and B — > XsC,'^i~ 
Monte Carlo event sample. 

The 18 modes selected in this analysis only capture about 53% of the full set of final 
states. Approximately 60% of the missing modes are due to final states with a K^j^ meson 
and their contribution can be determined from the modes. However, we need to account 
for the uncertainty in the fraction of modes with too many pions or kaons (two extra kaons 
may be produced via ss popping), as well as for modes with photons that do not originate 
from vr^ decays but rather from r], r^', etc. For final states with m{Xs) > 1.1 GeV/c^, we 
vary these fractions by ±5% per 7r°, ±20% for t], ±30% for Nt, > 5, and ±50% for rj' and 



Including systematic uncertainties, the measured branching fractions for m{i^i ) > 



where the first error is statistical and the second error is systematic. The combined B — >■ 
Xsi'^i~ branching fraction is the weighted average of the branching fractions for the electron 
and muon channels, where we assume the individual branching fractions to be equal for 
m(£+£^) > 0.2 GeV/c^. Table UITl summarizes the results of the analysis and lists both the 
statistical and systematic errors in the signal yields, the signal efficiencies and the branching 
fractions. 

The branching fractions for each Mx^ and bin are also measured, and summarized in 
Table HVl Figures a) and (b) show the distributions of the differential branching fractions 
as a function of (a) hadronic mass Mx^ and (b) = M|+^_ for electron and muon channels 
combined. 

SUMMARY 

Using a sample of 152 x 10^ T{AS) BB events, we measure the branching fraction 
for the rare decay B Xsi~^i~, where i = e or fi and Xs is a hadronic system is semi- 
inclusively reconstructed using 18 different hadronic states (with up to four pions). For 



others. 



0.2 GeV/c2 are 




(2) 
(3) 
(4) 
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TABLE III: Summary of results: signal yield (Ngig), statistical significance (Signif.), efficiency (e) 
and branching fraction (B). In the case of the signal yield and the branching fraction, the first 
error is statistical and the second error is systematic. In the case of the signal efficiency, the first 
error corresponds to uncertainties in detector modeling, BB counting, and Monte Carlo statistics, 
whereas the second error corresponds to the uncertainties in the signal model. 

Mode Nsig Signif. e (%) B (xlQ-^) 

XsC+e- 31.8 ±10.2 ±3.1 3.6 2.59 ± 0.18 1[J;3? 4.04 
36.3 ± 9.3 ±2.1 4.7 2.89 ± 0.14 l[J-j| 4.13 
Xsi+r 68.4±13.8±5.0 5.8 2.74 ± 0.13 tg-jo 4.11 



+ 1 so+"-*^" 
±0.83ir7'o 



TABLE IV: Branching fractions (B) for each bin of Mx^ and q'^. The first and second errors are 
statistical and systematic, respectively. 



bin 


B (xlO-7) 


bin 


B (xlO-7) 


MxA GeV/c2) 

[0.4,0.6] 

[0.6,0.8] 

[0.8, 1.0] 

[1.0,1.4] 

[1.4,2.0] 


3.75 ± 0.96l|];^^ 
0.36 ± 0.88+[!:j]^ 
6.65 ± 2.25l[!Ji 
10.50 ± 6.90t 
46.59 ± 23.37tlli\ 


GeV/c)2 
[0.04, 1.0] 
[1.0,6.0] 
[6.0, 14.4] 
[14.4,25.0] 


11.34 ±4.83t|^^ 
14.93 ± 5.04l|^^ 

7.32±6.14li;^^ 
4.18±1.17t°:^? 



m{i~^i ) > 0.2 GeV/c^, we observe a signal of 68.4 ± 13.8{stat) ± 5.0{syst) events and 
obtain a branching fraction of 

B{B Xstr) = (4.11 ± 0.83(stat) la7o(syst)) x 10"^ 

with a statistical significance of 5.8 a. 

TMs .suit . co,.,stent w.tK the .cent predict,™ by AU Q ou. p.ev.ous .c.us.ve 

B Xsi~^i~ measurement Q], and that of the BaBar collaboration within errors. 
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FIG. 5: Differential branching fraction as a function of (a) hadronic mass Mx^ and (b) = M'^ 



e+e- 



for electron and muon channels combined for data (points) and Monte Carlo signal (histogram). 
The outer (inner) error bars represent the total (statistical) errors. 
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